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ABSTRACT 

It is commonly assumed that ground-based gravitational wave (GW) instruments will not be sensitive to supermassive black 
holes (SMBHs) because the characteristic GW frequencies are far below the ~ 10-1000 Hz sensitivity bands of terrestrial 
detectors. Here, however, we explore the possibility of SMBH gravitational waves to leak to higher frequencies. In particular, if 
the high frequency spectral tail asymptotes to h(f) oc f~ a , where a < 2, then the spectral amplitude is a constant or increasing 
function of the mass M at a fixed frequency / ^> c 3 jGM. This will happen if the time domain waveform or its derivative exhibits 
a discontinuity. Ground based instruments could search for these universal spectral tails to detect or rule out such features 
irrespective of their origin. We identify the following processes which may generate high frequency signals: (i) gravitational 
bremsstrahlung of ultrarelativistic objects in the vicinity of a SMBH, (ii) ringdown modes excited by an external process that has 
a high frequency component or terminates abruptly, (iii) gravitational lensing echos and diffraction. We estimate the order of 
magnitude of the detection signal to noise ratio for each mechanism (i, ii, and iii) as a function of the waveform parameters. In 
particular for (iii), SMBHs produce GW echos of inspiraling stellar mass binaries in galactic nuclei with a delay of a few minutes 
to hours. The lensed primary signal and GW echo are both amplified if the binary is within a ~ 10deg(r/100M) -1 ' 2 cone behind 
the SMBH relative to the line of sight at distance r from the SMBH. For the rest of the binaries near SMBHs, the amplitude of 
the GW echo is ~ 0. l(r/100M) _1 of the primary signal on average. 



1. INTRODUCTION 

Gravitational wave (GW) observations are expected to open 
a new window on the universe within the decade. First 
generation GW detectors (LIGO, Virgo, TAMA, GEO) have 
been successfully commissioned, and the development of the 
next advanced-sensitivity ground based detectors (Advanced 
LIGO, Advanced Virgo, KAGRA) is well underway. How- 
ever, the construction of the planned space-based detector 
LISA/NGO may be delayed until after 2020. 1 

It is well known that the GW frequency 2 , generated by a 
purely gravitational process involving mass M, typically does 
not exceed M~ l . In particular, for quasicircular black hole 
(BH) inspirals, the spectrum peaks near twice the orbital fre- 
quency. As the separation shrinks, the GW frequency in- 
creases gradually to a maximum of order /j sco = (6 3//2 7rM) _1 = 
4.4Mg' mHz (where M 6 = M/(10 6 M Q )) corresponding to the 
innermost stable circular orbit. Then the BHs fall in quickly 
and the GW frequency reaches the fundamental quasinormal 
ringdown frequency / qnl - = (5.4irM)~ l = 12Mg'mHz. Simi- 
larly, hyperbolic, parabolic, eccentric, and zoom-whirl en- 
counters have a spectrum that peaks near f p ~ (27rr p /vp) _1 < 
lOMg'mHz, where r p > 3M and v p < 1 are the distance 
and relative speed at closest approach. These numbers cor- 
respond to nonspinning BHs, and vary by a factor of ~ 10 
for spinning binaries (Berti et al. 2009). This suggests that 
GWs are in the sensitive frequency band of terrestrial instru- 
ments (10-1000 Hz) only for low mass objects below M max ~ 
10 4 M Q (Flanagan & Hughes 1998; Amaro-Seoane & Santa- 
marfa2010). 

While it is clear that LISA would be very sensitive to GWs 
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1 http://lisa.nasa.gov, http://www.elisa-ngo.org 

2 We adopt geometrical units G = c = 1, and suppress the G/c 2 and G/c 3 
factors to convert mass to length and time units. 



emitted by SMBHs, we raise here the possibility whether 
ground based GW instruments could also detect massive 
sources. Since the instantaneous GW power output during 
a merger of compact objects is a constant multiple of c 5 /G ~ 
10 52 W, independent of M, which persists for a long time of 
order Af > 10M = 50Me s, it is not obvious a priori, whether 
a small portion of this enormous GW power leaking to high 
frequencies may be nonnegligible compared to the detector 
sensitivity level. 

We emphasize that the common assumption that black holes 
spacetimes cannot radiate at high frequencies / 3> M~ l is 
not true in general. In this paper, we review previous re- 
sults that show SMBHs are, at least in principle, capable 
of emitting or perturbing high frequency GW signals that 
LIGO/VIRGO may detect. We present general arguments 
that can be used to search for high frequency SMBH waves 
in the actual LIGO/VIRGO data, and motivate further theo- 
retical research to identify astrophysical mechanisms driving 
such modes. 

High frequency spectral tails, at frequencies well beyond 
the characteristic scale, M~ l , may be expected to follow cer- 
tain universal trends. A purely gravitational process involving 
BHs has only one dimensional scale: the mass, M. The mass 
scaling of the GW spectrum follows from dimensional analy- 
sis. The Fourier transform of the dimensionless GW strain for 
a source at distance D is most generally 

kf) = ^(Mf), (i) 

where ip is an arbitrary dimensionless function which may de- 
pend on dimensionless source parameters. 3 We are interested 
in whether there is any circumstance where the spectrum is 
an increasing function of M at fixed frequency / M~ l if 

3 To ensure a finite GW energy requires %l>(Mf) < (M/)~'/ 2 as / — > oo. 
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fixing all dimensionless parameters. Based on Eq. (1), this is 
possible if h(f) oc f~ a at high frequencies for a < 2. 

Similar effects, producing extended powerlaw spectral tails, 
are common in other fields of science and are often related 
to discontinuities in the signal or its derivatives. Examples 
are the "leakage" in data analysis when signals are abruptly 
terminated at the end of the sampling interval (Press et al. 
1992), or "termination waves" and "wind contamination" in 
geophysical gravity waves (Pulido & Caranti 2000). A sim- 
ilar 1 //" spectrum is also produced by stochastic processes 
in physics, technology, biology, economics, psychology, lan- 
guage and even music (for reviews, see Press 1978; Milotti 
2002). Here, we demonstrate that spectral GW tails associ- 
ated with (approximate) discontinuities have a universal high 
frequency tail, /"", where a is a positive integer, independent 
of the details of the waveform. Among these spectral profiles, 
Eq. (1) shows that an increasing function of mass requires 
a = 1 or 2. 

The distinct universal shape of these high frequency spec- 
tral tails allows to conduct searches for these GW burst 
sources irrespective of their origin. The relative arrival times 
and amplitudes of such signals in a network of detectors may 
be used to reduce false alarms and locate the source on the 
sky. 

If SMBHs were detectable, this would greatly increase 
the scientific reach of terrestrial GW instruments. Embed- 
ded in gas, some SMBHs shine as bright active galactic nu- 
clei (AGN), detectable with electromagnetic observatories. 
GW detections from SMBHs could serve to trigger searches 
for EM counterparts (Holz & Hughes 2005; Kocsis et al. 
2006), similar to counterpart searches for merging neutron 
stars (Nissanke et al. 2010, 2012; LIGO Scientific Collabo- 
ration et al. 2012; Abadie et al. 2012). An EM counterpart 
could be used to constrain cosmology (Schutz 1986; Dalai 
et al. 2006; MacLeod & Hogan 2008), extra dimensions (Def- 
fayet & Menou 2007) and the mass of the graviton (Kocsis 
et al. 2008). High frequency features, sensitive to small scale 
variability of the spacetime around SMBHs, may constrain 
the model of gravity. 

In this paper, we provide arguments which may guide fu- 
ture investigations to identify astrophysical mechanisms in- 
volving SMBHs that may produce high frequency GW spec- 
tra for LIGO/Virgo detections. First, we highlight the gen- 
eral connection between spectral tails and discontinuities us- 
ing Watson's lemma. Then, we discuss three examples where 
high frequency leakage may occur: ultrarelativistic encoun- 
ters, quasinormal ringdown modes, and GW echos. We high- 
light the need for further research in this direction to make 
more detailed assesments of these and other possibilities gen- 
erating high frequency spectral tails. 

2. THEOREM ON SPECTRAL TAILS 

Watson's lemma states that any piecewise infinitely dif- 
ferentiable, asymptotically exponentially decaying function, 
h(t), with a discontinuity in the rfi 1 time-derivative at an in- 
stant to, has an asymptotic spectral tail 4 \h(f)\ oc /~ (n+1) for 
/-»oo (Watson 1918). 

Watson's lemma can be generalized to smooth waveforms, 
which do not have an exact discontinuity in any derivative, but 

4 To avoid confusion, note that the high frequency GW tail discussed here 
is a spectral feature which has nothing to do with the late-time temporal GW 
tail of quasi-normal modes due to the scattering of GWs off the curved back- 
ground (Price 1972; Gleiser et al. 2008). 



for which the n th derivative undergoes a rapid transition within 
(comoving) time Afo <C M. Combining with the dimensional 
analysis argument, Eq. (1), 



M l ~ n 

«/>!<* ^ 



l+n ' 



for (Af (te r 1 »/»M ? - 1 , (2) 



where for a source at cosmological redshift z, luminosity dis- 
tance d h , M z = (l +z)M, At 0z = (1 +z)Af , and / = / em /(l +z) 
are the observed (cosmological-redshifted) quantities where 
/ em is the local emission frequency near the source. The di- 
mensionless constant of proportionality is independent of M, 
d h , and /. 

We conclude that the spectral amplitude may be an increas- 
ing function of the source mass M at high frequencies for fixed 
dimensionless parameters if n < 1, i.e., when the signal am- 
plitude or its derivative undergoes a rapid transition within a 
time Afo z <C 1 //. In the following we present specific exam- 
ples where this may be expected. 

3. GRAVITATIONAL BREMSSTRAHLUNG 

Gravitational bremsstrahlung during ultrarelativistic flybys 
of compact objects with a high Lorentz factor V is an exam- 
ple where the time-domain GW signal exhibits an approxi- 
mate discontinuity in the first derivative (Kovacs & Thorne 
1978). The previous dimensional arguments suggest that 
ground based instruments may be sensitive to GWs from such 
high speed flybys of objects near SMBHs. How large does the 
relative velocity have to be for detection? 

The GWs have been calculated for ultrarelativistic encoun- 
ters where b, the impact parameter, is sufficiently large to en- 
sure that the two objects move on nearly linear trajectories 
(Kovacs & Thorne 1978). The time domain GW signal is a 
"v-shaped" burst with a sharp discontinuity in the first deriva- 
tive and is beamed in the direction of motion. In the restframe 
of mass mi, 

Kf) « 4 ^ (1 + f )F [( 1 - v)K 2 (w) - ( 1 + v)K 2 (u)] 

in the forward direction at angles 9 <C T" 1 ' 2 , where T = 
(l-v 2 ) -1 / 2 is the Lorentz factor, v is the relative velocity. 
Here, ^(x) is a modified Bessel function, u = 2-Kfb/(vY), 
w = 2irfb/[v(l + v)T 2 ], q = mi/nix < 1 is the mass ratio, 
b is the impact parameter, and b = b/m\. Assuming that 
m\ ~ M ^> ni2, self-consistency requires b » 1 (see Death 
1978; D'eath & Payne 1992; Segalis & Ori 2001, and discus- 
sion below for the limit b ~ 0). We find that the signal extends 
into the Advanced LIGO/Virgo frequency band (/ > 10 Hz) if 

T > \%M l J 2 b x l 2 . Eq. (3) shows that the spectral amplitude 
is independent of M for fixed dimensionless parameters and 
fixed /, as it should, based on the dimensional arguments of 
Eq. (2). Figure 1 (bottom panel) shows that the amplitude 
may well exceed the Advanced LIGO/Virgo sensitivity level 
for SMBHs at 1 Gpc. 

Whether or not such ultrarelativistic objects exist near su- 
permassive black holes is unknown at present, but we argue 
that this possibility cannot be ruled out on theoretical grounds 
and may be tested observationally with Advanced LIGO and 
VIRGO. Ultrarelativistic emission with Lorentz factors be- 
tween r ~ 5 to 40 is observed in active galactic nuclei with 
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jets (Urry & Padovani 1995; Jorstad et al. 2001; Sikora et al. 
2005; Jorstad et al. 2005). These jets are believed to be 
launched by electromagnetic processes near a Kerr black hole 
(Blandford & Znajek 1977; McKinney & Blandford 2009). 
The total mass outflow rate during flares, inferred from obser- 
vations, reaches values comparable to the mass accretion rate 
of the BH, up to M Q /yr (Celotti & Ghisellini 2008). These 
observations suggest that at least energetically, the magnetic 
field and angular momentum of a Kerr SMBH is sufficiently 
large to accelerate a neutron star or white dwarf to ultrarela- 
tivistic velocities. Segalis & Ori (2001) considered the high 
frequency GW signal during the acceleration of a stellar mass 
localized blob of matter. In this case, h(t) ~ 4T 3 m/d h in 
the forward direction and 2Tm(l +cos6')/iiL at large angles 
9 ^> T" 1 . Here m is the mass of the accelerated object, and the 
frequency extends to / cut ~ V 2 / At, where Af is the accelera- 
tion timescale in the SMBH restframe. 

Secondly, a merging SMBH binary may gravitationally ac- 
celerate stellar mass compact objects to ultrarelativistic ve- 
locities, if they happen to be in the vicinity of the merging 
binary (van Meter et al. 2010). The latter possibility is sup- 
ported by several arguments. First, a widely separated SMBH 
binary gravitationally traps stellar objects in mean motion res- 
onances from the stellar cluster surrounding the binary. These 
objects are analogues of Trojan satellites in the Solar System 
corotating with Jupiter. An accretion disk around the binary 
may also supply objects to mean motion resonances by cap- 
turing stars crossing the disk or forming stars in the outskirts 
of the disk and transporting them radially inwards close to 
the binary (Subr & Karas 1999; Karas & Subr 2001; Ayliffe 
et al. 2012). As the SMBH binary inspirals, objects locked in 
a mean-motion resonance are dragged inwards to small radii 
down to 10 Schwarzschild radii (Schnittman 2010; Seto & 
Muto 2010, 2011). Numerical relativity simulations by van 
Meter et al. (2010) have shown that test particles surround- 
ing a SMBH binary may be accelerated to high Lorentz fac- 
tors during the final merger of the binary. If so, the GW 
bremsstrahlung of these objects due to the SMBH may be can- 
didates for LIGO/VIRGO detections. Future studies should 
investigate this process in more detail for a variety of initial 
conditions, SMBH binary mass ratios and spins, and stellar 
mass gravitating particles to determine the likelihood of such 
detections. 

4. BLACK HOLE RINGDOWN 

A perturbed BH behaves as a damped harmonic oscillator, 
and emits GWs (for a review, see Berti et al. 2009) 



h + (t) 
h x (t) 



k 



-t/n 



cos(2irf k t) 
sin(27r/*?) 



(4) 



where + and x denote the polarization, t > is time from the 
end of the external perturbation, ft is the quasinormal ring- 
down frequency, 1 /Vfc is the decay rate, and Si^ are spin-2 
weighted spheroidal harmonics which set the angular distri- 
bution of the emission with root mean square (Att)~ x I 2 . Here k 
is a discrete index which runs over all of the ringdown modes 
(n,£,m) analogous to the principal, angular momentum, and 
magnetic quantum numbers of the hydrogen atom. While /j. 
and T/t are universal properties of BHs which depend only 
on their mass and spin, the Ak amplitude coefficients, called 
quasinormal excitation factors, are sensitive to the details of 
the perturbation. The excitations factors can be computed as 
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FIG. 1 . — Characteristic spectral amplitude for two waveform families. Top: 
Quasinormal ringdown at distance 100 Mpc for mass M = 1O 6 M0, amplitude 
A = O.lr, and decay time t/M = 10 (left) and 0.001 (right curves), respec- 
tively. Red and blue curves correspond to different conventions at t < 
(see § 4). Red and blue curves represent driving processes that terminate 
abruptly and gradually, respectively. Bottom: Ultrarelativistic flybys with 
r = 50, impact parameter b = 20M, and mass ratio q — 10~ 5 for M — 10 6 
(red) and 10 5 Mq (blue) at 1 Gpc. Solid green curves show the sensitivity 
of LISA (left) and Advanced LIGO (right). Advanced LIGO may detect 
weakly damped SMBH ringdown modes excited by a process that terminates 
abruptly and ultrarelativistic encounters. 

a convolution of an external source with the Green's function 
of the SMBH (Leaver 1986; Andersson 1995; Dorband et al. 
2006; Berti & Cardoso 2006; Baker et al. 2007; Berti et al. 
2010; Dolan & Ottewill 201 1; Hadar et al. 201 1). 

To compare with the sensitivity levels of GW instruments, 
one needs to calculate h(f), the Fourier transform of the wave- 
form. For a damped harmonic oscillator, the high frequency 
tail depends on the start of the waveform. Two conventions 
are used in the literature discussing the detectability of ring- 
down waveforms (Echeverria 1989; Finn 1992; Flanagan & 
Hughes 1998; Berti et al. 2007, 2006). In the Echeverria-Finn 
(EF) convention, the signal is assumed to vanish for t < 
and follow Eq. (4) for t > 0. In the Flanagan-Hughes (FH) 
convention, to avoid an unphysical discontinuity at t = 0, the 
signal is extrapolated symmetrically for t < and the overall 
amplitude is renormalized by 2" 1 / 2 to keep the total GW en- 
ergy fixed. However, the signal used in the FH convention is 
also discontinuous in the first derivative, which has an impact 
at high frequencies. EF and FH introduced these conven- 
tions for convenience to calculate the detectability near the 
characteristic ringdown frequency of the signal, but here we 
explore the high frequency asymptotics of these conventions 
and identify the corresponding physical conditions. Phys- 
ically, an EF spectrum is appropriate in the case where the 
QNM is excited by an instantaneous pulse. The FH spectrum 
represents excitation processes that build up the mode con- 
tinuously, but generate a sudden change in the derivative of 
the waveform as the excitation is terminated abruptly at t = 0, 
letting the BH ring freely thereafter. For waveforms excited 
smoothly on timescales rjt, we introduce a third convention 
which is continuous in all derivatives, by replacing e~'l Tk by 
V2sech(t /r k ) = V8/(e' /Ti + e-'/ Ti ) in Eq. (4). 
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The Fourier transform of the ringdown waveform can be 
obtained analytically in the three conventions for each mode 
k comprising the waveform, independently. In the limit of 
asymptotically large frequencies / » max(/t, -r^T 1 ), we get 
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where (2t = 7r7>/fc is the quality factor of the damped oscil- 
lator. Equations (5-7) show that the asymptotic spectral tail 
may enter the detector band even if / > M" 1 . For fixed /, the 
spectral GW amplitude scales very differently with M. Since 
Ak, Tk and 1 are all proportional to M for fixed k, we find 
that \hH<j)\ ocM, \hf x (f)\ and \h^(f)\ are independent of 
M, and ocM" 1 , while \hf* h (f)\ and |/z£ ch (/)l are asymp- 
totically exponentially suppressed for / >• M~ x . In the later 
case, the peak spectral amplitude is proportional to Aujk oc M 2 
for a fixed k at frequencies / <~ r^ 1 . At / <~ t^T 1 for some A:, 
the orientation-averaged one-sided dimensionless GW spec- 
tral amplitude per logarithmic frequency bin for this signal is 

{2fhf ch } ~ (4/5) 1 / 2 (47r)- 1 / 2 (2- 1 / 2 7 r) A/^l- These mass scal- 
ings also follow from dimensional analysis once the spectral 
shape is specified (§ 1) and are consistent with the general 
expectations from Watson's lemma (§ 2). The EF and FH 
spectra become invalid beyond a frequency / > Af" 1 , where 
At represents the characteristic timescale over which the exci- 
tation occurs generating the approximate discontinuity in the 
signal amplitude or derivative, respectively. 

The top panel of Figure 1 shows the angular averaged spec- 
tral amplitude forM = 10 6 M Q , d L = 100 Mpc, and A k = O.lr. 
Different line types correspond to different conventions at the 
start of the waveform as labelled. Two sets of curves peaked at 
<~ 0.01 and 30Hz represent r = 10M (weakly damped modes) 
and 10~ 3 M (highly damped modes), respectively. The GW 
signal leaks into the LIGO/Virgo frequency band for FH- 
type excitation processes both for weakly and highly damped 
modes. If the excitation occurs smoothly over timescale r 
(blue curves), then the signal may enter the LIGO/Virgo band 
if r = 1O 3 M0, but not if r is much larger. 5 Comparing to 
the minimum frequency and sensitivity level of Advanced 
LIGO/VIRGO, we conclude that SMBH ringdown waveforms 
may be detectable on average 6 to <~ 100 Mpc if the quasinor- 
mal modes are excited to an amplitude exceeding Ak > 0.01 r. 

Going beyond the statement that some QNMs may have a 
potentially detectable high frequency tail, it would be desir- 
able to identify conceivable astrophysical processes exciting 
these modes to detectable levels. The Ak excitation factors of 
QNMs are explored only for a very limited number of config- 
urations to date (Leaver 1986; Andersson 1995; Dorband et al. 
2006; Berti & Cardoso 2006; Baker et al. 2007; Berti et al. 
2010; Dolan & Ottewill 201 1; Hadar et al. 201 1). For this 
discussion, we distinguish three classes of ringdown modes. 

5 Note that the X polarization of the waveform is suppressed at high fre- 
quencies because the signal vanishes at t = in this case, see Eq. (4). Only 
the + polarization is detectable at high frequencies for smooth excitation pro- 
cesses for highly damped modes (blue dashed curve). 

6 averaging over polarization and direction 



• High-^ modes. For of a Schwarzschild BH, f n g m w 
m c /(2Tr) and 1 /r nem w nfl c , where fl c = 27 _1 /2 M _1 is 
the orbital frequency at the light ring (Press 1971). 7 
These modes may be identified with energy leakage 
from the photon sphere (Goebel 1972; Ferrari & Mash- 
hoon 1984; Dolan & Ottewill 2011; Yang et al. 2012). 
They ring in the LIGO/VIRGO band if 10 Hz < f nlm < 
1000 Hz. 

• Weakly damped modes. Consider modes with f n ^ m < 
M~ l , 1 /T„£ m <~ M~ l . Although, the ringdown frequency 
is well outside the LIGO/VIRGO frequency range in 
this case for large M, a rapidly changing external driv- 
ing process on timescale Af < 0.1 s may lead to high 
frequency tails for these modes (see e.g. EF and FH 
above). This is analogous to driving a harmonic oscil- 
lator with a process that has a component well above 
the resonant frequency. Further studies are necessary 
to see if the high frequency spectral amplitude of the 
QNM may be sufficiently large for any specific driving 
process in practice. 

• Highly damped modes. The decay rate of very high 
overtone modes may be arbitrarily large, lim„ 1 / T n i m = 
oo for all I and m, which implies a high frequency spec- 
tral tail. Although the quality factor of these modes 
approaches zero, their excitation may still be sufficient 
to reach levels required for detection as stated above 
(A k > 0.01t at lOOMpc). For a damped oscillator 
driven by a stationary external process of amplitude less 
than unity, Ak < Tk/(Mfk). These modes are detectable 
to 100 Mpc with Advanced LIGO if/* < 100 M~ l and 
r <~ 0.1 s (see the right blue curve on the top panel of 
Figure 1 and discussion above), which is satisfied for 
all modes with I < 3000. 

We highlight two astrophysical processes which may gen- 
erate high frequency QNMs. First, a merging SMBH binary 
gravitationally accelerates particles and nearby stellar mass 
compact objects to ultrarelativistic velocities which may then 
escape or merge with the SMBH (van Meter et al. 2010, and 
see § 3 above). Ultrarelativistic collisions with a SMBH 
(Davis et al. 1972; Berti et al. 2010) excite high-^ ringdown 
modes. The radiated energy scales as Eg oc l/£~ dE , where 
dg ~ 2 for radial infall and 1 for grazing captures, implying 
that h(f) - / _3/2 to f~ 2 for high /. Thus, similar to ul- 
trarelativistic gravitional bremsstrahlung, the high frequency 
spectral GW amplitude of ultrarelativistic capture sources is a 
constant or increasing function of mass for a fixed mass ratio 
(see Eq. 1). 

Secondly, the high frequency gravitational effects of a stel- 
lar (or intermediate) mass (Mb) black hole binary may excite 
high frequency QNMs of the SMBH. Stellar mass compact 
binaries form during close encounters and merge due to GW 
emission (O'Leary et al. 2009; Kocsis & Levin 2012), or they 
may be driven to merger due to the Kozai effect (Wen 2003; 
Antonini & Perets 2012; Naoz et al. 2012). The QNM exci- 
tation function by these sources has a high frequency com- 
ponent and the driving force terminates on a short timescale 
proportional to M\,. Based on the arguments outlined above, 
it is then plausible to expect a high frequency ringing of the 

7 Similar expressions hold for Kerr BHs, see Ferrari & Mashhoon (1984); 
Yang et al. (2012); Keshet & Ben-Meir (2012). 
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SMBH. Furthermore, the high frequency GWs emitted by 
such stellar (or intermediate) mass compact binaries are scat- 
tered by the SMBH, generating GW echos. The scattered 
waves represent high-i? ringdown modes, which we discuss 
in § 5 below. Further studies are necessary to examine the ex- 
citation factors and the detectability of the signal to assess the 
likelihood of detecting these waveforms. 

5. GW LENSING ECHOS AND DIFFRACTION 

Consider the GWs emitted by an inspiraling stellar mass 
BH binary source in a galactic nucleus near a SMBH. High 
frequency GWs are scattered by the gravitational field of the 
SMBH (e.g. Futterman et al. 1988) producing GW echos 
(Dolan & Ottewill 2011; Zenginoglu & Galley 2012). The 
GW frequency is not effected to leading order in the WKB 
limit if the frequency of the incoming wave exceeds M" 1 . The 
scattered GWs, or GW echos, are completely analogous to 
the relativistic images of electromagnetic waves for rays go- 
ing around the SMBH in opposite directions and/or multiple 
times. However, GW echos may be easier to distinguish 
from astrophysical foregrounds than the relativistic images of 
EM sources since high frequency GW sources are sufficiently 
rare, have a short duration, and a distinct frequency evolution 
(O'Leary et al. 2009; Kocsis & Levin 2012). Furthermore, 
GW detectors are sensitive to the signal phase and amplitude, 
i.e. square-root of the intensity, which is suppressed much 
less than the EM intensity (see Fig. 2 below). 

The GW echo arrives at the detector after the primary signal 
with a delay corresponding to the projected light-travel time 
between the object and the source, 

Ar = (l-cosa)r- 14h x (1 -cosa)M 6 (r/10 4 M), (8) 

where r is the distance between the binary and the SMBH and 
a is the deflection angle. Note, that the probability distribu- 
tion of inspiraling stellar mass BH binaries in galactic nuclei 
formed by GW emission during close encounters increases to- 
wards the SMBH (O'Leary et al. 2009). A signficant fraction 
(~ 30%) of sources are within r < 10 4 M for which the GW 
echos occur within hours of the primary signal. 

We calculate the amplitude of the GW echo for a source 
at radius r from the SMBH in the geometrical optics limit 
using the Virbhadra & Ellis (2000) lens equation valid for ar- 
bitrary deflection angles. We direct the reader to that paper 
for details (see also Bozza 2008). For a sanity check, we also 
calculate the scattering cross-section of plane waves imping- 
ing on the SMBH, which approaches the lensing cross-section 
in the large deflection angle limit if r ^> M. The amplitude 
of the GW echo may be larger or smaller than the primary 
GW signal of the source depending on the relative position 
of the source, the SMBH, and the line-of-sight as shown by 
Figure 2. The lensing effect may strongly increase the inten- 
sity of the wave if the lens and the source are nearly colinear 
along the line of sight (Einstein 1936). For r <C D, the de- 
flection angle corresponding to the Einstein radius for strong 
magnification is a E ~ (4r/M) _1 / 2 , i.e. (11°, 3.6°, 1.1°) for 
r = (10 2 , 10 3 , 10 4 )M, respectively. For an isotropic distribu- 
tion of sources, this implies that (2,0.2,0.02)% of sources are 
within the Einstein radius for the three values of r, respec- 
tively. For larger deflection angles, the GW amplitude de- 
creases as h tx ho/(ra 2 ), where ho is the GW amplitude with- 
out lensing. The peak at a ~ tt (parallel backscattering) is 
known as the glory effect of geometrical optics (Futterman 
et al. 1988). The figure shows that the amplitude of relativis- 
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FIG. 2. — Gravitational lensing magnification of a high frequency GW 
source by a SMBH as a function of deflection angle a. Colored curves use 
the thin lens equation, black dotted lines show the plane wave limit. The three 
sets of curves correspond to different source distances r from the SMBH as 
labeled, different line styles represent relativistic images of different order. 
Solid curves show rhe primary signal; dashed and dotted-dashed curves show 
the amplitude of the first GW echo. The peak at small deflection angles cor- 
responds to strong lensing at the Einstein radius ole = (4r/M)~'/ 2 , the second 
peak at 180° represents the glory phenomenon in geometrical optics. Addi- 
tional peaks appearing at a = rnr, where n is an integer, are not shown beyond 
n = 1. The GW frequency is assumed to greatly exceed M~ l . Diffraction 
fringe patterns (not shown) appear with spacing 8a ~ (34Mf) and sharper 
features are smoothed. 

tic images may be significant even for large deflection angles 
if r < 10 3 M and the primary signal has an SNR larger than 
~ 50. 

Figure 2 assumes / ^> M~ l and neglects diffraction effects 
due to the interference of rays going around the SMBH in 
opposite directions and the QNM of the SMBH (Matzner 
et al. 1985; Nakamura 1998). This is a good approximation 
for typical large deflection angles of stellar mass GW inspi- 
rals at the typical distances r ^> 100M, since the timescale 
for the increase of the GW frequency is much less than AT 
(Eq. 8). However, interference effects may become signifi- 
cant for sources that happen to be behind or in front of the 
SMBH close to the line of sight. This smears out the mag- 
nification peak and produces a GW fringe pattern in am- 
plitude and phase on the angular scale 5a ~ l/(34/M) ~ 
5.9 x 10" 4 x M^C/yiOHz)" 1 for a nonspinning SMBH (sim- 
ilar expressions with different prefactors hold for spinning 
SMBHs, see Dolan 2008). 

The frequency and direction dependence of the GW diffrac- 
tion patterns introduce interesting time variations in the sig- 
nal. The diffraction pattern scale shrinks in time for inspiral- 
ing binary sources as / increases (Dolan 2008). Furthermore, 
the relative motion of the lens and the source lead to a vary- 
ing GW amplitude analogous to the microlensing events of 
electromagnetic signals (Mao & Paczynski 1991) and phase 
shift. Assuming that the source moves along a circular or- 
bit around the SMBH with a Keplerian angular velocity, f^K, 
it travels across the diffraction fringe within 5tf = Sa/Q.^ ~ 
9.3s x (r/10 3 M) 3 / 2 (//100Hz)" 1 . This is much less than the 
source lifetime for inspiraling stellar mass binaries. Based on 
the strong magnification rates estimated above, we conclude 
that a small fraction of sources (< 2%) might exhibit these 
features. 

If measured, the GW echo may be used to learn about 
the parameters of the SMBH beyond its distance from the 
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merging stellar-mass binary. Time dependent effects of GW 
diffraction patterns of an inspiral signal and/or the microlens- 
ing effect discussed above carry information about the mass 
of and spin of the SMBH. The polarization of the scattered 
waves relative to that of the primary signal is also sensitive to 
the SMBH spin (Dolan 2008). 

6. DISCUSSION 

Ground based instruments may, in principle, be sensitive to 
GWs emitted by SMBHs with arbitrarily large M at frequen- 
cies / » M~ l . Dimensional analysis shows that the spectrum 
h(f) is a growing function of mass or mass-independent at 
large frequencies if h(f) oc or f~ 2 , respectively. Watson's 
lemma shows that such a spectrum corresponds to a wave- 
form having an abrupt change in amplitude or in its first time- 
derivative within a short time (Af < 1 //). Discontinuous GW 
shock fronts with arbitrarily small Af are self-consistent, ex- 
act solutions of general relativity (Vaidya 1951; Aichelburg & 
Sexl 1971; and see also cosmic strings for which h(f) oc /" 1//3 , 
Damour & Vilenkin 2000, 2005), so one should not discard 
this possibility a priori. 

Ground-based GW detectors could provide observational 
constraints on the event rates of such sources independent of 
their origin. Unlike other GW burst sources, these sources 
have a well-defined broadband spectral shape, h{f) oc f~ l or 
f~ 2 , allowing to develop optimal search algorithms for their 
efficient detection. These signals may be distinguished from 
instrumental glitches if it shows up coincidentally in all of the 
different instruments in a detector network (i.e. the two LIGO 
sites, Virgo, and Kagra) with a consistent spectral amplitude. 

We have discussed three hypothetical examples involving 
SMBHs potentially capable of generating high-frequency GW 
signals: ultrarelativistic sources, QNM ringdown waveforms 
excited externally by high frequency processes or those 
that terminate abruptly, and the scattering and diffraction of 
GWs. First, an ultrarelativistic stellar mass object mov- 
ing in the gravitational field of a SMBH emits gravitational 
bremsstrahlung. The signal derivative undergoes a rapid tran- 
sition on an observed timescale f~^ t ~ Af /T 2 , where T is the 
Lorentz factor and Af is the acceleration timescale propor- 
tional to the impact parameter, b. Then for M~ l -C / < / cut , 
the GW spectrum h(f) oc f~ 2 is proportional to the mass 
of the object but independent of the mass of the SMBH, 
M, extending well into the sensitive range of ground based 
GW instruments if T > 5.7 M 5 b. Here M 5 = M/10 5 M© and 
b = b/M. Second, SMBHs, if perturbed, ring as damped 
harmonic oscillators, and emit ringdown GWs. We have 
shown that weakly damped ringdown waveforms leak into 
the LIGO frequency band if the excitation process has a high 
frequency component or terminates within 0.1 s. These GWs 
may be detectable for arbitrary M 3> 100M Q , with Advanced 
LIGO/Virgo-type instruments to lOOMpc if the GW ampli- 
tude satisfies Ak > O.Olr, where t is the decay time of the 
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QNM. However, detailed calcalations of excitation factors 
are necessary to determine whether or not this level may be 
reached for specific realistic configurations. 

Finally, lensing and diffraction effects of high frequency 
GWs by SMBHs may also generate detectable features. In 
particular, we expect a 0.02%-2% of LIGO/Virgo sources in 
galactic nuclei (e.g. inspiral and merger of stellar mass black 
hole or neutron star binaries) could be strongly magnified by 
the central SMBH. Diffraction limits the maximum magni- 
fication of these sources. The rest of the high frequency 
sources in the vicinity of SMBHs produce weak GW echos 
due to large deflection angle scattering by the SMBH. The 
GW echo amplitude is typically between 0.1-10% of the pri- 
mary signal for these sources. The time delay between the 
primary and scattered signals in this case is typically tens of 
minutes. 

It remains to be seen whether these processes occur in real- 
ity. To date, numerical simulations of binary mergers did not 
find evidence for waveforms having such sharp features for 
quasicircular mergers (Baker et al. 2007), except for spurious 
GWs related to the improper choice of initial conditions (Field 
et al. 2010). However these simulations are limited by reso- 
lution to comparable mass binaries (see Lousto & Zlochower 
201 1 for the most unequal mass ratio) and cannot resolve high 
order ringdown modes. Recently, Berti et al. (2010) carried 
out simulations for head-on, grazing, and off-axis collisions 
with nonzero initial velocity. They found that the cutoff fre- 
quency is shifted to larger values for higher Lorentz factors 
and multipole moments. Future studies should investigate 
other possible examples, involving SMBHs. These may in- 
clude high frequecy ringdown resonances in extreme mass ra- 
tio inspirals or three-body encounters. High frequency GWs 
of SMBHs, may also arise in alternative theories of gravity 
or cosmology due to the self-gravity of GWs (Fodor & Racz 
2003; Kaloper 2005). 

One should also be aware of these spectral effects when 
constructing time domain templates matching the inspiral, 
merger, and ringdown waveforms (Babak et al. 2007; San- 
tamarfa et. al. 2010). A high-frequency powerlaw spectral 
tail is sensitive to a discontinuity in the derivative of the 
time-domain waveform at the matching boundaries. We have 
shown that the high frequency QNMs are very sensitive to 
the initial conditions for both the EF and the FH conventions. 
This theorem may also be useful for debugging numerical 
codes. 
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